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Abstract. Imbibition phenomena have been widely used experimentally and theoretically to study the ki- 
netic roughening of interfaces. We critically discuss the existing experiments and some associated theoretical 
approaches on the scaling properties of the imbibition front, with particular attention to the conservation 
law associated to the fluid, to problems arising from the actual structure of the embedding medium, and 
to external influences such as evaporation and gravity. Our main conclusion is that the scaling of moving 
interfaces includes many crossover phenomena, with competition between the average capillary pressure 
gradient and its fluctuations setting the maximal lengthscale for roughening. We discuss the physics of 
both pinned and moving interfaces and the ability of the existing models to account for their properties. 

PACS. 47.55. Mh Flows through porous media - 05.40+j Fluctuation phenomena, random processes, and 
Brownian motion - 68.35.Ct Interface structure and roughness 
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1 Introduction 

A considerable amount of effort has been spent in the 
field of kinetic roughening over the last two decades. Apart 
from technological interests in crystal growth, a good part 
of the fascination for this field comes from the possibility 
of describing many different types of interfaces by a few 
distinct universality classes, in terms of scaling exponents 
and scaling functions 

However, if the theoretical and numerical aspects of 
the field are extremely rich and varied, the experimental 
backing of these ideas is quite lacking. This is for instance 
the case in kinetic roughening in crystal growth due to 
the great variety of atomistic processes. Much of the ex- 
perimental attention has rather concentrated on simpler 
one-dimensional "toy-systems" , in particular driven inter- 
faces in random media. A definite advantage of these sys- 
tems is that the interface configuration is directly observ- 
able in the experiment (while, it must often be deduced, 
from some probe-surface interaction, in crystal growth). 
Even these simple systems are challenging. They are most 
commonly grouped in universality classes described by the 
Kardar-Parisi-Zhang g (KPZ) or Edwards- Wilkinson [§ 
(EW) equations, with either quenched or thermal noise, 
depending on the driving regime. In the case of quenched 
noise the exact low-dimensional scaling behaviour is how- 
ever still somewhat contradictory j|,^,^,||j|]). 

Examples of experimental studies in these systems are 
slow combustion fronts jlO| (for which KPZ scaling was re- 
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cently demonstrated), shock fronts, fluid-gas interfaces in 
Hele-Shaw cells |l]jl|Jl|,|l|,^ , or in paper [|6l0Jl|,|l|,|ol||l 
and fracture surfaces |25| such as one-dimensional fracture 
lines 1^,^ , to name a few. In this paper, we concentrate 
on a particular system, namely the spontaneous imbibition 
of a porous medium by a liquid |l|,y,|l8|jl||^,|l|,|||2|,|| . 
Our goal is to critically review the experiments and theo- 
ries that exist in the literature and to indicate a direction 
for further investigation and comprehensive understand- 
ing of imbibition front roughening. The key questions are 
(i) whether imbibition should present any kind of scaling 
behaviour at all, and (ii), if so, under which conditions 
can one of the - possibly several - scaling regimes become 
observable. 

Our motivations in writing this paper are twofold. First, 
the fluid dynamics aspect of imbibition is itself very com- 
plex, and should not be ignored in any discussion of the 
statistical fluctuations of the interface. We feel that this 
has not been properly done so far. Secondly, we present in 
a companion paper the results of a line of investigation of 
spontaneous imbibition based on a phase field formalism 
[ P8[ ; the present paper is intented to lay down the aspects 
which we believe essentials to any models of imbibition. 

We start in Section ^ by reviewing the macroscopic 
properties of imbibition and show that even the simple 
propagation of an imbibition front may have several dy- 
namical regimes, depending on the design of the exper- 
iment. We also show that the macroscopic progression 
of the interface will have a very strong influence on the 
roughening process. Section || outlines the models that 
have been proposed for imbibition experiments as well as 
their predictions. In Section ^, we discuss the existing ex- 



2 



M. Dube, M. Rost and M. Alava: Conserved Dynamics in Imbibition 



perimental data, and point along the way to several fea- 
tures missed by these experiments. Finally we conclude 
with a discussion and proposals concerning future experi- 
ments. 



2 Macroscopic Features of Imbibition 

Although extremely familiar to researchers working in the 
field of flow in porous media, the details of imbibition 
would seem to be relatively unknown to the statistical 
physics community. It is generally defined in reference to 
two-fluid flow in porous media, and corresponds to the 
displacement of the lesser wetting fluid by the more wet- 
ting one [ p9||30[ |. Notice that this definition is irrespective 
of whether the flow of the fluids is spontaneous or induced 
(eg., by a pump). In this work, we will restrict ourselves 
to the case of spontaneous imbibition. The flow of the flu- 
ids is thus driven solely by capillary forces, with gravity 
and/or evaporation being the only external influences on 
the fluids' motion. 



2.1 Capillary Rise 

The simplest example of spontaneous imbibition is capil- 
lary rise: a part of a capillary tube, of radius R, is im- 
mersed into a reservoir, exposed to an ambient atmo- 
spheric pressure Pq. We assume that the fluid wets the 
capillary so that a meniscus, described by the surface 
z — h{x, y,t) = is formed. At equilibrium, this surface is 
characterised by a contact angle 9, obtained from Young's 
law |l| 

Jig cos 6 = Jsg -Ish (1) 

where jij is the surface tension between the phases i and j, 
and s, / and g stand for solid, liquid and gas respectively. If 
the meniscus is in motion, the contact angle differs from its 
equilibrium value [3^ , but we neglect this effect here. The 
important point is that the curvature causes a pressure 
difference AP = Pc ^ Pn ~ P{z = /i") = 2-figCOse/R 
where 7 is the surface tension of the liquid-gas interface 
and R/ cos 9 is the curvature of the meniscus. We refer to 
Pc as the capillary pressure. 

The motion of the fluid, of density p and viscosity 77, 
can be treated within the assumption of Poiseuille flow 
p3t , i.e. the full Navier-Stokes equation is replaced by the 
simpler Stokes equation for an incompressible fluid under 
the influence of the gravitational force pg. As usual, the 
pressure field is found from Laplace's equation V^P = 0, 
but with boundary condition P{z — 0) — Pq and P{z — 
h) = Pq — Pc- (We can neglect corrections to the pressure 
field close to the meniscus, if its height is a lot larger than 
the radius of the tube). This yields P{x,y,z) = P{z) = 
Pq — Pcz/h{t), and it is straightforward to obtain the pro- 
gression of the interface, 

dh K / heq \ 

dt = [W) ) ^ ^ 



a classical result derived by Washburn Q and Rideal . 
Washburn's and Rideal's equation includes the permeabil- 
ity of the tube k = R'^/8 and the equilibrium height of the 
meniscus heq = Pc/ pg- Studies of capillary rise including 
the inertial term of the Navier-Stokes equation show 
that Eq. (^) is essentially correct, with the exception of 
very short times. Notice also that this equation of motion 
also neglects completely the problem of the actual motion 
of the contact line between the gas-solid and liquid phases 




The (transcendental) solution of Eq. (g) has the fol- 
lowing asymptotic properties: Defining /i(to) as the initial 
height of the column and r^q = h^qri / K{pg)'^ as an equi- 
libration time, for low heights h ^ heq (and also in the 
absence of gravity where h^q = 00) the rise is of the form 

h\t)-h\to) = — (t-to). (3) 
V 

The equilibrium height is approached by h exponentially 

h{t) ^ /ie,(f-e-*/^-). (4) 

This was examined experimentally by Washburn [ p4| 
and Rideal |^ as well as numerous others, who found very 
good agreement between theory and experiment. Capil- 
lary rise in a liquid-liquid system was also examined by 
Mumley et al. |37| . They confirmed Washburn's result in 
the case of perfect wetting, but reported discrepancies for 
systems with non-zero contact angle. In this case, the rise 
was slower that t^/^, a behaviour attributed to the motion 
of the contact line itself Q . 

2.2 Capillary Rise in Porous Media 

With the basic capillary rise phenomenon understood, we 
can now turn to spontaneous imbibition in porous media. 
The flow of liquids in porous media is generally described 
in terms of Darcy's equation p9| , |30{ 

(Q) = -pA^(VP-pg), (5) 

where (Q) is the average mass of fluid transported per unit 
time through the cross-section A of the porous medium, 
and K and 77 are the average permeability and viscosity 
respectively. Darcy's equation arises from an averaging 
procedure of the porous medium and ignores all details 
on length scales smaller than the pores. If we assume the 
porous medium to be homogeneous, such that the per- 
meability is a constant kq independent of the fluid con- 
centration, we can still solve Laplace's equation for the 
pressure and obtain a coarse-grained pressure gradient, 
VP = zPc/h{t), with Pc an effective capillary pressure, 
and h the average height of the fluid column. This is of 
course valid only when the notion of an interface is it- 
self well defined, i.e. it should not be too "fuzzy" up to 
macroscopic scales. Under these assumptions, and with 
the identification (Q) = pAdh/dt, Darcy's equation leads 
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directly to the Washburn-Rideal result, Eq. (g) with a 
dynamical rise given by Eqs. (^) and (||). 

This form of Darcy's equation has been used to study 
many experiments of fluid propagation in porous media, 
including some fibrous materials (see below). The agree- 
ment is no t a lways perfect - we shall come back to this 
in Section 2.1. At this point, a more complete system of 
equations [30|3^ would only obscure the physics. 



2.3 Evaporation 

Considering the typical experimental setups for imbibition 
front propagation ||2^, evaporation effects should also be 
included to Washburn's equation. As far as we are aware, 
no detailed studies of fluid motion through thin porous 
media including the effects of evaporation have been per- 
formed. A natural assumption is to introduce an evapo- 
ration rate proportional to the area of the fluid exposed 
to air, i.e. the net loss of fluid mass per unit time due 
to evaporation is Qe = —2£ pLh{t) where L is the lateral 
width, and 5 is a phenomenological evaporation rate per 
unit area. Thus Eq. (0) is modified to 



dh 
'dt 



1 (I 

n \ 



'Pc_ 

rj \ h 



P9 



eh 



(6) 



where e — 2£p is the evaporation rate. This form obviously 
neglects changes in the concentration profile towards the 
surface, as well as the evaporation at the interface itself, 
which is both justified for thin media such as a sheet of 
paper. 

An immediate consequence of Eq. (^ is an equilibrium 
height depending on the evaporation rate 



heq 



2h„ 




(7) 



where hg = Pc/ pg and h(, = (Pci^j tr\)^l'^ are respectively 
the equilibrium heights in the cases where only gravity or 
evaporation have a significant infiuence, with a crossover 
defined by 2hg ~ /ig. When gravity can be neglected, the 
height h(t) behaves as 

h^{t) = h^(ta) e~2-(*^*«) ^hl (l - e-2^(*"*°)) (8) 

Well below h^. the rise still follows h(t) ^ t^^^, and it again 
approaches exponentially the equilibrium height. 



2.4 Addition of a Solution 

Many experiments on the roughness of the interface in an 
imbibition context have studied the behaviour of a dye, 
added to the pure fluid ]l6| , p0p^ . In the literature on 
porous media this is referred to as the hydrodynamic dis- 
persion phenomenon [|30|. At the simplest level, it can be 



treated by the introduction of a dye concentration field 
K(x,t) advected by the macroscopic velocity field v. 



OK 
'dt 



V • VK = DV^K, 



(9) 



where £) is a diffusion constant (in some cases, it may be 
necessary to introduce a diffusion tensor Dij ||3^ ) . 

Two remarks on Eq. (^ are in order. First, in every- 
day life, it is easily seen that the fluid front is always faster 
than the dye front. This clogging phenomenon, which may 
be due to a smaller permeability for the dye, may easily 
be incorporated phenomenologically by introducing a con- 
stant A < I in the front of the convective term, i.e. 



V • Av ■ VK. 



(10) 



The second point concerns the stopped interface. If the 
pinning of the interface, at a distance h^, from the reser- 
voir, is due to evaporation, there will nevertheless be a 
fluid motion, with approximate velocity Ve ~ nPc/Tjhf,, in 
order to compensate the losses. This implies that more 
and more dye particles will be brought to the fluid in- 
terface, thus creating a band of high dye concentration, 
analogous to the coffee rings examined recently |^^. The 
increase of the width of this band with time can be used as 
an alternative way of measuring the relevant macroscopic 
parameters. 



2.5 Validity of Macroscopic Description 

We now discuss the validity of the macroscopic results for 
the specific case of spontaneous imbibition in paper. Ordi- 
nary paper is made out of wood fibers and, in many cases, 
chemical additives and filler materials like talc and clay, 
arranged in a disordered structure. Not only is there a 
wide distribution of pore sizes, with a high effective tortu- 
osity but the surface of the network is extremely uneven. 
This gives rise to all the standard problems in defining 
a static/dynamical contact angle for inhomogeneous sub- 
strates. 

Darcy's equation, Eq. (^, is relatively well established 
(to the point that it is even often referred to as a law) for 
general porous media. However, paper, as well as other 
fibrous materials present the peculiarity that the fiber 
structure may be modified by the contact with the liq- 
uid, a phenomenon known as swelling |^^. Cellulose fibers 
show a great affinity to water and can absorb large quan- 
tities in millisecond timescales, giving rise to concomitant 
changes in fiber volume and pore structure (this is how- 
ever not the case for many organic fluids and oils). Thus 
one's intuitive picture of fibers as capillary tubes is false: 
the pore structure is highly non-trivial and in some cases 
time-dependent. There are two different effects that play 
a role: the volume to be filled with liquid increases, and 
the flow resistance of the pore network changes. 

There can also be an exchange of liquid between the 
inside of a flber and the "surface" pores, which compli- 
cates enormously the fluid flow, since there are no well 
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defined "structures" (either the pores, or the fibers) re- 
sponsible for the capillary forces |^ 42|. This is of course 
a very serious impediment to any kind of attempt at a 
"microscopic description" of the structure, say in terms of 
a percolation network p3[ . 

A more serious problem with water is the fact that it 
is not actually known whether a pressure balance argu- 
ment is always valid. On short time scales (smaller than 
a few seconds), it has been proposed that front penetra- 
tion would proceed in pores first with the establishment 
of a prewetting layer through diffusion Bq] , in which case 
the front position could advance as i with k varying 
up to unity. Experiments have demonstrated that clear, 
Washburn-like behavior can be obtained in conditions that 
amount to basically no interaction (swelling, prewetting) 
between the material and the penetrating liquid . Note 
that minute applications of chemicals, e.g. during paper 
manufacture, may induce drastic changes in the effective 
viscosity or surface tension of the invading liquid. 

There are however some advantages to paper. The very 
high permeability of the fibers will reduce the formation of 
overhangs in the interface. Darcy's equation cannot treat 
trapped air bubbles in bulk porous media |38 , which in 
paper only play a minor role: paper is thin and the pores 
should be connected well to the surface. Moreover spon- 
taneous imbibition is slow and allows more time for the 
removal of overhangs and trapped bubbles. 

In conclusion, one would believe that a well defined 
fluid-air interface should exist, and that Darcy's equation 
should be valid on length scales larger than the interface 
width, provided that interactions between the liquid and 
the fibers are minimal. In that sense, a full hydrodynam- 
ical treatment of the problem may not be necessary (as 
noted in [2^ ). In many cases, it may also be necessary 
to introduce a concentration dependent permeability, i.e. 
K = k{p) [Q. In any case a Washburn approach may 
be a first step. In that respect, recent experimental work 
with paper and deionised water, showing front propaga- 
tion consistent with Washburn's equation is quite encour- 
aging 0. 



2.6 Statistical Fluctuations of the Interface 



So far, we have considered a flat interface in a completely 
homogeneous medium. In general, the disorder in the pa- 
per causes fluctuations which will accumulate to roughen 
the advancing front. In the standard picture of kinetic 
roughening [|l]j|] the fluctuations of the interface are cor- 
related up to a distance $|[(t), a lateral correlation length, 

increasing in time as ^|| ~ t^^^ and described by the dy- 
namical exponent z. The vertical extent of the interface 
fluctuations, its width W, is related to ^|| through the 
roughness exponent W ^ Thus the width increases 

with time a,s W ^ t'^ with /3 — x/z, until the interface 
fluctuations have saturated, ie., until £,\\{t) = L, which 
defines the crossover time ty L^. 



The initial increase of the width and eventual satura- 
tion are comprised in a (Family- Vicsek) scaling form 



W{L,t)^L^f{t/L^) 



(11) 



with the scaling function f{u) ~ for u 1 and f{u) ~ 
const for 1. 

The spatial height difference correlation function (of 
the qth moment) 



(12) 



Gg{r,t)^{\h{r + r',t)-h{r',t)n 

generally scales similarly to the total width, i.e. Gq ^ 
for small r, and approaching a constant as r — > oo. How- 
ever, care has to be taken: In the case of anomalous scal- 
ing the height differences for fixed |r| do not saturate 
with time, and G„ shows the local roughness exponent 



Xioc 48|. The height difference distribution for fixed |r 
may also have a long tail, causing intermittent or "turbu- 
lent" jumps in the height configuration in which case each 
moment has a different exponent Xq, ie., the interface ex- 
hibits multiscaling |Q. Also a temporal height difference 
correlation function 

Cg{t) = {\h[v,t+t')-h{v,t')-{h{t+t')~h{t'))\'i)l''i (13) 

can be considered (e.g. in [^ ), which is of particular use in 
the case the system exhibits time-translational invariance. 

To our knowledge none of the works studying interface 
fiuctuations in imbibition have checked for anomalous (as 
determined by the higher moments of the correlation func- 
tions) types of scaling behaviour |l|,|l|jl|Jl|,|^,|l|j2|,||,|4 
although this is in principle difficult, since an exact diag- 
nosis for higher moments g > 2 is severely hampered by 
large statistical fluctuations. These concepts are never- 
theless important when the randomness of the medium is 
described with random-field disorder |p8| . 

Another peculiarity of imbibition is that, apart from 
the standard kinetic lateral correlation length £,\\{t), aris- 
ing from the accumulating history of fluctuations, another 
lateral length scale is to be seen: local conservation of the 
fluid determines a lateral length related to the average 
motion of the interface. Intuitively, it is clear that such 
a length scale must exist, since fluctuations ahead and 
behind the average position of the imbibition front have 
respectively a slower and faster instantaneous velocity. To 
make a quantitative argument, we introduce an effective 
surface tension 7*, representing the energy cost of a curved 
air-liquid interface on a macroscopic scale [js^ . Any cur- 
vature at the interface - the typical size of fluctuations 
being W vertically and ^x laterally - modifles the pres- 
sure by an amount AP — 7*VF/^x (the Laplace pressure 
effect). This should be compared with the difference in 
the pressure field across the same vertical distance W due 
to the pressure gradient derived in Section p.l| , namely 
AP ~ PcW/H, from which we obtain 



(14) 



relating the parallel length scale ^x to the height of the 
interface. Beyond this length scale we do not expect any 
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correlated roughness of the interface, since those fluctua- 
tions would be suppressed by the overall gradient in the 
pressure field Pc/H. 

Let us recall once again that the interface continuously 
slows down p^ , without gravity nor evaporation H ^ 
t^/^. Thus the slow increase of ~ H^^^ ~ i^^"* with 
time leads to an increase of the width W ^ with f3 — 
x/4. However, we stress that this increase is conceptually 
different from the increase of the lateral correlation length 
with a dynamic exponent z in standard models of kinetic 
roughening. 



3 Models of Imbibition 

In parallel with experimental work, many theoretical mod- 
els of imbibition have been developed. These models fall ei- 
ther in discrete (cellular automata) or continuous classes. 
Most of the theoretical discussion on imbibition was done 
in terms of the DPD model, which was shown to belong 
to the same universality class as the quenched KPZ equa- 
tion. Other types of models or continuum equations have 
also been proposed. 



average capillary pressure. In that sense, any overhangs 
will eventually be removed, although the question of the 
different time scales involved in the dynamics is certainly 
extremely complex [ p8| . 

The DPD model in its simplest form then probably 
describes well the statistical properties of the pinned in- 
terface. It is however difficult to believe that this model 
can describe the whole dynamical motion of the interface, 
since it does not account for liquid conservation. It already 
fails by predicting a constant average interface velocity. As 
we saw in section [2], one of the principal features of im- 
bibition is a constant slowing down of the front, a result 
intimately related to fluid conservation. It can be argued 
that the DPD refers speciflcally to the motion of the dye, 
but the rise of the dye front is necessarily bounded by the 
liquid-air interface, and the conservation law governing the 
motion of the fluid must be reflected on the motion of the 
dye particles. From that point of view, simply replacing 
V in Eq. (^ certainly would be incorrect. 

It is interesting to show how the global behaviour of the 
modifled DPD model may be related to the macroscopic 
Washburn behaviour, Eq. (||). The density of open cells 
1 — p{h) can be thought of as an effective force acting on 
the interface, similar to the gradient in pressure on fluid 
motion. Since the Washburn Eq. is by essence a dissipative 
equation, we can associate, away from pinning 



(15) 



3.1 Cellular Automata Models 

The earliest theoretical investigations of imbibition phe- 
nomena were based on the Directed Percolation Depinning 
(DPD) model (l|,|l^,|§. It is a cellular automaton, with 
space discretized in cells which are either blocked or wet- 
table, with the blocked cells being a fraction p. The par- 
ticularity of the model is that overhangs in the interface 
are removed as soon as they occur, which thereby intro- 
duces anisotropy in the interface. Growth takes place by 
invasion of available cells until the interface comes across 
a percolating directed path of blocked cells. Clearly, the 
statistical properties of the interface are related to the 
properties of the percolating path. Close to the percolation 
threshold pc — 0.47, the directed percolating path is char- 
acterized by the parallel and perpendicular length scales 
■ClI ~ b~Pc|~''" and ^ Ip^PcI^"^ , with values ~ 1.7 
and i^±_ ~ 1.09. so the width W{L) - - L^^/^w , yield- 
ing X — 0.63. 

To treat evaporation, the DPD model was modified by 
the introduction of height dependent fraction of blocked 
cells, i.e. p = p{h) ~ h/h^ with an equilibrium height 
Hq ^ (Vp)~^ defined by p(/io) = Pc In this model, 
the interface is described by the standard DPD roughness 
exponent up to a length Zx ^ [\/p)~'^/^'^''° after which it 
saturates to a constant value Wsat (Vp)"'''. The identi- 
fication Wsat ^± then yields 7 = J^±/(l + J^j_) ^ 0.52, a 
value comparable to the experimental results. 

As already noted ||2l[, there is a physical justification 
to the main feature of the DPD model, namely the erosion 
of overhangs. It is plausible that propagation of the fluid 
in a direction parallel to the reservoir, if it has a chance to 
occur, will be favoured over the motion of the fluid away 
from it, since this does not result in any change in the pinning transition, as well as a change of morphology of 



The assumption of constant gradient thus effectively cor- 
responds to an evaporation rate. It is also possible, at the 
macroscopic level, to include gravity, and to mimic the 
slowing down of the interface. Note however that we pre- 
dict a pinning height he ~ e^^^'^ while the assumption of 
Amaral et al. implies ^ (Vp)^^, since they neglect the 
capillary driving term. 

Another model of evaporation was introduced by Ku- 
mar and Jana j2^. The model is essentially similar to 
the DPD model, with the modification that each cell is 
not necessarily "full" or "empty" , but may contain many 
"subunits" . Evaporation is modeled by a loss of n subunits 
in the transfer between cells. It is found that below a criti- 
cal Uci , the interface propagates indefinitely. Between this 
value and a second critical loss rate ?ic2, they observe an 
interface behaviour similar to the one observed by Ama- 
ral et al. 1 20 , but with value x = 0-5 and 7 = 3.0. For 
n > nc2, this scaling regime breaks down and the rough- 
ness exponent becomes n-dependent. These results were 
backed experimentally, but not in any consistent way. 

It is rather difficult to believe that indefinite front 
propagation under evaporation (obtained for n < Ud) is 
physically realistic, nor seems the way of including evapo- 
ration convincing. In this model, the concentration of fluid 
molecules decreases continuously from a maximum Nq to 
a value of at the interface, which seems physically unre- 
alistic. 

We terminate by presenting invasion models based on 
the random field Ising model which predict a de- 
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the interface at some length scale related to the capillary 
pressure. It seems however unsuited for spontaneous im- 
bibition, because any advance of the front requires the 
increase of an externally applied pressure. 

A variation of this idea, introduced by Sneppen |50| , 
allows invasion always at sites of lowest resistance. An 
interesting aspect of this model is temporal multiscaling 
PH , due to avalanche motion 1 52 5^ . Avalanches are also 
present in the process of spontaneous imbibition, although 
the consevation law imp oses a natural cutoff on their size 
and distribution |]l5| , ^8[ . It is nevertheless interesting that 
a similar lack of temporal scaling is also seen in the phase 
field model of imbibition, thus surviving the introduction 
of a conservation law. 



4 Experimental Efforts in Imbibition 

Most of the experiments concerned with front propagation 
in random media seem to produce a self-afRne interface, 
but the numerical value of the exponents often bears lit- 
tle resemblance to the standard universality classes and/or 
associated models. Imbibition of porous media would seem 
like an obvious first experimental choice, since the asso- 
ciated time and spatial scales are easily accessible in the 
laboratory. Many such imbibition experiments were per- 
formed, but it turns out that there is very little in common 
between the different experiments. This is partly due to 
the fact that different experiments had different goals, but 
it is also a reflection of the great complexity of the pro- 
cesses involved in imbibition. 



3.2 Continuum Description 

The DPD model (in its original version without "evap- 
oration") belongs to the same universality class as the 
quenched KPZ equation (see e.g. Chapter 10 in ^ and 
references therein, as well as [^) 

= i.V'hix,t)-^iVh{x,t)f + ^{h{x,t),x). (16) 

A related equation with multiplicative noise has been in- 
troduced by Csahok et al. |^ in order to model the random 
porosity of the medium. Since these equations are local at 
the interface, they do not contain a conservation law for 
the fluid. 

With these considerations in mind, a continuum equa- 
tion was introduced by Zik et al. following Krug and 
Meakin |54]. In Fourier space, it has the form 



dt 



Kp. 



|fc|-^ 
V hoit) 



-Vkit) 



(17) 



where k is the wave- vector, hnCt) is the mean height of the 
interface at time t (cf., Eq. (g|)), and rjk is the noise term, 
assumed to be annealed, with correlations {i]k{t)i]k' {t')) ^ 
6k+k'S(t—t'). This equation is essentially similar to the one 
derived by Krug and Meakin for the roughness of stable 
Laplacian fronts with noise and corresponds to the leading 
term of a Saffmann- Taylor analysis of the problem. Note 
that |fc| is nonlocal in space. 

The other difficulty with respect to spontaneous im- 
bibition, the quenched nature of the disorder, remains. 
As long as the interface sweeps fast enough through the 
medium, the disorder acts as annealed, time dependent 
noise , although the ques tion of the effective noise cor- 
relator is far from trivial |28| — Eq. |l^, as it stands, gives 
the unphysical results of an interface that never saturates. 
In any cases, at late times under slow average interface 
propagation in a continuum model one is confronted with 
both the difficulties of a nonlocal interface equation with 
quenched disorder. 



4.1 Experiments on Pinned Interfaces 

The earliest statistical ph ysic s imbibition experiments were 
done by Buldyrev et al. |l|,|l^,0 and Family et al. 
The first experiment was performed with a dye solution in 
a vertical capillary rise setup. The rising front moved from 
the reservoir and the roughness of the interface developed 
during the process. Eventually, the dye front stopped, due 
to gravity and/or evaporation (no dynamical measure- 
ments, either micro- or macroscopic, were done in this ex- 
periment), and the roughness of the pinned interface was 
measured. The main experimental conclusion of the paper 
was a roughness exponent with value x = 0.63, consistent 
with the DPD model (see above, Section H). However, the 
length scale over which the scaling behaviour was observed 
is extremely small. For a sheet of paper of total lateral ex- 
tent L — 40 cm, C{1) ~ /-^ only for length scales I smaller 
than Imax ~ 1 cm, after which it levels off to a constant 
or to a logarithmic function of I |5^. Since this scaling 
region is only a few times larger that the length scale of 
of the fibers themselves, it is not at all clear whether the 
exponent x really is a universal value, or simply results 
from the microscopic fiber structure. As a comparison, 
the scaling region for the paper burning experiment [ |lO| 
was from 1.5 cm to about 10 cm, for a total length 30 cm 
[ p6[ . A similar experiment was done in a three dimensional 
sponge-like material [ p^ . Again, the stopped interface was 
observed, yielding a roughness exponent x^^''^ ~ 0.5, a re- 
sult consistent with the higher-dimensional version of the 
DPD model. No other experimental details were however 
given. 

In a further set of experiments, Amaral et al. stud- 
ied the role of evaporation in more detail. They again 
considered the stopped interface, but changed the height 
where pinning occured through the evaporation rate (pre- 
sumably by modifying the humidity during the experi- 
ment, although this is not specified). Their main result is 
that the width of the pinned interface is related to the 
pinning height Hp (and thus to the evaporation strength) 
through a new exponent 7, namely w^at ~ H^. The ex- 
periments gave a value 7 = 0.49, which - as reported in 
Section ^ - can also be obtained from a modified version 
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of the DPD model [gOj . (We will come back to this in Sec- 
tion O). On the other hand, Kumar and Juma |22 also 



performed an experiment in presence of various evapora- 
tion rates and claimed that the roughness exponent was 
not universal but stronly dependent on the evaporation 
rate, i.e. x — x{^)^ although no systematic experimental 
investigation of this breakdown of universality was done. 



4.2 Experiments on Moving Interfaces 

The experiment of Family et al. was performed in an hor- 
izontal capillary setup with water only, and the position 
of the air-water interface was recorded both temporally 
and spatially p8| . The main experimental results were an 
average interface progression h ^ , with 6 ~ 0.7, and 
a self-afHne interface described by a Family- Vicsek scal- 
ing relation, and characterized by the exponents /3 = 0.38 
and X ~ 0.76. The first result h{t) t'^-'' is certainly in- 
triguing. However, it should be noticed that a reservoir 
was placed underneath the piece of paper, in order to pre- 
vent evaporation. It is of course possible that this caused 
condensation instead, thus increasing the velocity of the 
front. It is also in line with the results of non- Washburn 
water penetration in paper |Q . Also here the spatial scal- 
ing regime was rather small, for distances below Imax ~ 2 
cm for a 40 cm wide sheet of paper. It should also be 
noticed that these experiments were made with Chinese 
paper, having a thickness of only a few fiber layers, which 
may also influence the results. 

The temporal scaling of the interface was studied in de- 
tails by Horvath and Stanley . Moving a piece of paper 
such that the interface always remained at a fixed distance 
H above a reservoir, they found a power law behaviour for 
the time correlation function C2[t) ~ with /3 = 0.56. 
Another result is that the velocity V at which the the piece 
of paper must be moved towards the reservoir in order to 
keep the interface at a constant H varies as, V ^ H~^^, 
where Q — 1.6. Notice how this implies an interface prop- 
agating such that H{t) ^ t^/^, since V = dH/dt. This is 
slower than 5—1/2 expected from Darcy's equation, but 
consistent with our earlier discussion. Unfortunately, the 
scaling behaviour of the interface as a function of space 
was not discussed and/or measured in this otherwise very 
careful work. 

It is certainly interesting that the value of [3 is con- 
stant for all heights considered. The different heights did 
however affect the time at which the time correlation func- 
tions saturated, and the values at which it did so. Horvath 
and Stanley found for C{t) a scaling form 



(18) 



where f{y) is a scaling function such that f{y) ^ for 
y ^ 1 and f{y) ^ const for j/ 3> 1. The values of the 
independent exponents were 0l = 0.48 and Or = 0.37 

Another experiment was performed by Kwon et al. 
psf . They deposited a paper towel on an inchned glass 
plate and followed the capillary rise of a dye solution, giv- 
ing a mean interface dynamical propagation h(t) i'^ '^^. 



Again, two scaling regimes were present in the saturated 
width; on small length scales (< 2 cm), x — 0-67 while 
on larger length scales (up to 20 cm) x ~ 0.2. Within the 
simple Family- Vicsek picture they obtained f3 = 0.24 on 
the short lengthscale regime. 

Finally, Zik et al. p4| performed an horizontal capil- 
lary front experiment. They obtained rough interfaces only 
with highly anisotropic paper, for which they found x = 
0.4. For isotropic paper, the roughness was at best loga- 
rithmic. It is remarkable, that the scaling for the anisotropic 
paper was observed through a large range of length scales, 
not only up to fiber length. 



4.3 Analysis of the Experiments : Influence of Fluid 
Conservation 

As the previous paragraphs show, there is a large amount 
of experimental work on the subject of spontaneous imbi- 
bition and front roughening, but a coherent picture does 
not emerge naturally. The results are quite contradictory 
and difficult to analyse properly. 

We have already mentioned a few experimental diffi- 
culties inherent to imbibition in paper, the most impor- 
tant being the complicated flow proflle of a fluid inside the 
paper, due to fiber swelling and dissolution. Certainly all 
experiments suffer from these difficulties at least to some 
degree. Apart from the work of Family et al. ||l^ , all ex- 
periments are consistent with an average front propagat- 
ing roughly according to the Washburn-Rideal result, Eq. 
(||). They don't agree perfectly, but definitely reflect the 
influence of the conservation of fluid on the propagation 
process. 

One consistent feature to emerge from the experiments 
is an interface which has developed roughness only on very 
short length scales. A partial remedy to this problem is, 
as done by Amaral et al. |^^, to introduce a length scale 
Ip above which the interface saturates. It their experiment 
and model (see also Section ||) it becomes visible in the 
width the interface takes at the height where it gets pinned 
by evaporation, and they establish the relation w ^ 
ip. 

The idea of such a length scale has of course a wider 
range of applicability than in the context of imbibition 
with evaporation. Indeed, the experiment of Horvath and 
Stanley , done with little influence of evaporation, and 
no discernable influence of gravity, showed the existence of 
a similar length scale, in this case, related to the velocity 
V at which paper was being pulled down (or, equivalently, 
to the average height of the interface). 

We have already shown that such a length scale (H), 
connected to the existence of a conservation law existed. 
It can provide a natural connection between the work of 
Amaral et al. and Horvath and Stanley, provided that 
one is willing to revise the role of evaporation. If the 
sole role of evaporation is to stop the macroscopic pro- 
gression of the interface, with little or no influence on 
the statistical fluctuations of the interface itself, then the 

width w{t) - C^(-ffp) H}^^ and therefore 7 = x/2. 
This length should also influence the setup examined by 
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Horvath and Stanley p8[ . Indeed, it provides an alterna- 
tive explanation to the velocity (or alternativley height) 
dependence of the saturated value of the time correlation 
function, since 

Clir^^)- ^ (19) 

which yields C2(t ^ oo) - _ y-'^^^. Note that 

as long as < L, the correlation function C2{t) is in- 
dependent of the total length of the system. Since only 
one length of paper was considered experimentally, it is 
unfortunately impossible to check this hypothesis. 

However, in order to demonstrate convincingly the ex- 
istence of such a length scale (or of any other length scales) 
and the scaling below it, it is imperative to extend the 
region of scaling. We fail to see how a sub-centimeter de- 
scription of the roughness might be considered universal. 
It is quite possible that the crossover fx associated with 
paper really is in the sub-centimeter range, in which case, 
the alternative way to check for scaling is in the scaling of 
the width. Even in the freely rising case, the total width 
of the interface should show some early time power law 
behaviour, with some exponent (3. 

One should also note the fact the structure of paper 
does have short-scale power-law correlations These 
are manifest in the two-point-correlation function of the 
areal mass density of the paper, with the cut-off of the 
effective power-law correlations extending in practice to a 
few times the fiber length. It is by no means clear what 
role such correlations play in imbibition, since penetra- 
tion of the liquid may take place both inside the sheet, 
and on the outer surface. Also, it is not clear how these 
correlations map into correlations in the hydraulic con- 
ductivity or local permeability. They however remind of 
the fact that scaling exponents established on lengthscales 
that are also relevant from the microstructural viewpoint 
should be taken with a grain of salt. Still, it is quite possi- 
ble that the effective noise correlator at the interface will 
display power law behaviour, as in the case of forced flow 
through Hele-Shaw cells pof . 



5 Summary and Perspective 

In this paper, we have tried to extract a coherent picture of 
all previous studies about interface roughening in sponta- 
neous imbibition experiments. Concerning theoretical un- 
derstanding, the DPD picture plausibly reproduces the 
properties of the pinned interface. We strongly feel that 
not enough attention has been paid to the macroscopic 
dynamical features and their connection to the underly- 
ing microscopic structure. The flow of fluid in paper is 
influenced by many microscopic mechanisms, and a study 
of the average front velocity already indicates which ones 
may or may not be relevant. The main experimental char- 
acteristic of imbibition certainly is the slowing down of 
the front, a phenomenon not seen in other systems. Some 



very simple predictions concerning evaporation and grav- 
ity would already allow a solid comparison between the 
expected behaviour and the experimental results. 

An important reason to consider the fluctuations of the 
propagating front is the presence of a lateral length scale 
fx , for which experimental evidence already exists | pO| , pl| , 
and which has a very strong influence on the roughening 
process. 

Finally, we believe that a complete description of im- 
bibition (including the dynamic interface) must include 
a conservation law for the whole fluid, not simply at the 
interface level. This is for two reasons, first the average in- 
terface ought to slow down continuously, and secondly the 
dynamics of the interface need to have long range inter- 
actions. Also, any model must deal the quenched nature 
of the noise in a suitable way. These requirements make 
it impossible to describe the interface roughening with a 
local equation of motion. To answer these criteria, we in- 
troduce in the next paper a phase-field model of im- 
bibition, which incorporates explicitely a conservation law 
for the liquid in a disordered medium, and thus produces 
the correct macroscopic physics. 

We finish this overview by pointing out a few exper- 
imental directions. The use of "ordinary" paper for the 
medium creates some problems that are usually neglected. 
In this respect it comes to mind that other liquids than wa- 
ter or water/ink mixtures might well prove advantageous, 
in that they are less susceptible to evaporation and/or do 
not interact with the fibers themselves. Another option is 
to monitor the roughness of the front using an Hele-Shaw 
cell. This eliminates the problem of fiber-liquid interac- 
tion. 



6 Acknowledgements 

We wish to acknowledge interesting and stimulating dis- 
cussions with J. Kertesz, T. Ala-Nissila, K. R. Elder, S. 
Majaniemi, and J. Lohi. This work has been supported by 
the MATRA program of the Academy of Finland. 



References 

1. See eg., A.-L. Barabasi and H.E. Stanley, Fractal Concepts 
in Surface Growth, (Cambridge University Press, Cam- 
bridge, 1995), 

2. J. Krug, Adv. Phys. 46, 139 (1997). 

3. M. Kardar, G. Parisi and Y.-C. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 

4. S.F. Edwards and D.R. Wilkinson, Proc. Royal Soc. Lon- 
don A 281, 17 (1982). 

5. H. Leschhorn, Phys. Rev. E 54, 1313 (1996). 

6. T. Nattermann, S. Stepanow, L.-H. Tang, and H. 
Leschhorn, J. Physique II 2, 1483 (1994). 

7. H. Leschhorn, T. Nattermann, S. Stepanow and L.-H. 
Tang, Ann Physik 6, 1 (1997). 

8. O. Narayan and D.S. Fisher, Phys. Rev. B 48, 7030 (1993). 

9. Z. Csahok, K. Honda, E. Somfai, E. Vicsek, and T. Vicsek, 
Physica A 200, 136 (1993). 



M. Dube, M. Rost and M. Alava: Conserved Dynamics in Imbibition 



9 



10. J. Maunuksela, M. Myllys, O.-P. Khknen, J. Timonen, N. 
Provatas, M.J. Alava, and T. Ala-Nissila, Phys. Rev. Lett. 
79, 1515 (1997). 

11. M.A. Rubio, Phys. Rev. Lett. 63, 1695 (1989). 

12. V.K. Horvath, F. Family, and T. Vicsek, J. Phys. A 24, 
L25 (1991). 

13. S. He, G.L.M.K.S. Kahanda, and P.-z. Wong, Phys. Rev. 
Lett. 69, 3731 (1992). 

14. T. Delker, D.B. Pengra, and P.-z. Wong, Phys. Rev. Lett. 
76, 2902 (1996). 

15. A. Dougherty and N. Carle, Phys. Rev. E 58, 2889 (1998). 

16. S.V. Buldyrev, A.-L. Barabasi, F. Caserta, S. Havlin, H.E. 
Stanley, and T. Vicsek, Phys. Rev. A 45, R8313 (1992). 

17. A.-L. Barabasi, S.V. Buldyrev, S. Havlin, G. Huber, H.E. 
Stanley, and T. Vicsek, in Surface Disordering : Growth, 
Roughening and Phase Transitions, R. JuUien, J. Kertesz, 
P. Meakin and D.E. Wolf (Nova Science, Commack, 1992). 

18. F. Family, K.C.B. Chan, J. G. Amar, in Surface Disorder- 
ing : Growth, Roughening and Phase Transitions, R. Jul- 
lien, J. Kertesz, P. Meakin and D.E. Wolf (Nova Science, 
Commack, 1992). 

19. S.V. Buldyrev, A.-L. Barabasi, S. Havlin, J. Kertesz, H.E. 
Stanley, and H.S. Xenias, Physica 191 A, 220 (1992). 

20. L.A.N. Amaral, A.-L. Barabasi, S.V. Buldyrev, S. Havhn, 
and H.E. Stanley, Phys. Rev. Lett. 72, 641 (1994). 

21. V.K. Horvath and H.E. Stanley, Phys. Rev. E 52, 5166 
(1995). 

22. P.B.S. Kumar and D. Jana , Physica A 224, 199 (1996), 
see also cond-mat/9509047, 

23. T.H. Kwon, A.E. Hopkins, and S.E. O'DonneU, Phys. Rev. 
E 54, 685 (1996). 

24. O. Zik, E. Moses, Z. Olami, and L Webman, Europhys. 
Lett. 38, 509 (1997). 

25. E. Bouchaud, J. Phys. Cond. Matt. 9, 4319 (1997). 

26. J. Kertesz, V.K. Horvath, and F. Weber, Fractals 1, 67 
(1993). 

27. T. Eng0y, K.J. Maal0y, A. Hansen, and S. Roux, Phys. 
Rev. Lett. 73, 834 (1994). 

28. M. Dube, K.R. Elder, M. Rost, M.J. Alava, S. Majaniemi, 
and T. Ala-Nissila, to appear Eur. Phys. J. B (2000). See 
also M. Dube et al., Phys. Rev. Lett. 83, 1628 (1999). 

29. A.E. Scheidegger, The Physics of Flow through Porous Me- 
dia, MacMillan Co, New- York (1957). 

30. M. Sahimi, Rev. Mod. Phys. 65, 1393 (1993). 

31. J. S. Rowlinson and B. Widom, Molecular Theory of Cap- 
illarity, Clarendon Press, Oxford (1982). 

32. P.G. de Gennes, Rev. Mod. Phys. 57, 827 (1985), other 
references may be found in, eg., J.F. Joanny and L. Leger, 
Rep. Prog. Phys. 55, 431 (1992). 

33. L.D. Landau, Fluid Mechanics (Pergamon Press, NY, 
1979). 

34. E.W. Washburn, Phys. Rev. 17, 273 (1921). 

35. E.K. Rideal, Philos. Mag. 44, 1152 (1992). 

36. J. Szekely, A.W. Neumann, and Y.K. Chuang, J. Colloid . 
Interface Sci. 35, 273 (1971). 

37. T.E. Mumley, C.J. Radke, and M.C. Wilhams, J. Colloid 
Interface Sci. 109, 398 (1986), and ibid 109, 413 (1986). 

38. R. Hilfer, Phys. Rev. E 58, 2090 (1998). 

39. R.D. Deegan et al.. Nature 389, 827 (1997). 

40. J. A. Bristow, Svensk Papperstidning 20, 645 (1971). 

41. J.S. Aspler, S. Davis, and M.B. Lyne, J. Pulp and Paper 
Sci. 13, J55 (1987). 



42. I. Pezron, G. Bourgain. and D. Quere, J. Colloid Interface 
Sci. 173, 319 (1995). 

43. R. Lenormand, E. Touboul, and C. Zarcone, J. Fluid Mech. 
189, 185 (1988). 

44. T. Gillespie, J. CoU. Science 13, 32 (1958), and J. Coll. 
Science 14, 123 (1959). 

45. See, eg., D.E. Eklund, and P.J. Salminen, Tappi Journal 
69 (9), 116 (1986). 

46. O. Zik, T. Kustanovich, E. Moses, and Z. Olami, Phys. 
Rev. E 58, 689 (1998). 

47. J. Krug, Phys. Rev. Lett. 72, 2907 (1994). 

48. J.M. Lopez, M.A. Rodriguez, and R. Cuerno, Phys. Rev. 
E 56, 3993 (1997). 

49. N. Martys et al., Phys. Rev. B 44, 12294 (1991). 

50. K. Sneppen, Phys. Rev. Lett. 69, 3539 (1992). 

51. K. Sneppen and M.H. Jensen, Phys. Rev. Lett. 71, 101 

(1993) . 

52. Z. Olami, I. Procaccia, and R. Zeitak, Phys. Rev. E 49, 
1232 (1994). 

53. H. Leschhorn, and L.-H. Tang, Phys. Rev. E 49, 1238 

(1994) . 

54. J. Krug and P. Meakin, Phys. Rev. Lett. 66, 703, (1991). 

55. This can be seen easily from Fig. (2) of Ref. jl^ ]. 

56. On length scales smaller than 1.5 cm, there may be signs of 
scaling, with a x ~ 0.7, but again, it is difficult to attach 
any universal significance to length scales of the order of 
the fiber length. 

57. Horvath ands Stanley also included the total system size L 
in a scaling hypothesis, but since only one paper size was 
used, their hypothesis could not be fully tested. 

58. As far as we are aware, only the experimental groups of 
Refs. [^, and have realised this physical fact, but 
the implications of it were not discussed in any elaborate 
ways. 

59. N. Provatas, M.J. Alava, and T. Ala-Nissila, Phys. Rev. 
E54, R336 (1996). 

60. V.K. Horvath, F. Family, and T. Vicsek, Phys. Rev. Lett. 
67, 3207 (1991). 



